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ABSTRACT 

The nature and properties of dark matter are both outstanding issues in physics. The universal 
character of gravity implies that sclf-gravitating compact configurations of dark matter are not only 
possible, but likely to be spread throughout the universe. The astrophysical signature of these objects 
may be used to probe fundamental particle physics, or even to provide an alternative description 
of compact objects in active galactic nuclei. Here we discuss the most promising dissection tool of 
these configurations: the inspiral of a compact stellar-size object and consequent gravitational-wave 
emission. The inward motion of this "test probe" encodes unique information about the nature of the 
central, supermassive dark matter configuration. When the probe travels through some compact dark 
matter profile we show that, within a Newtonian approximation, the quasi- adiabatic evolution of the 
inspiral is mainly driven by dark matter accretion into the small compact object and by dynamical 
friction, rather than by gravitational-wave radiation-reaction. These effects circularize the orbits and 
leave a peculiar imprint on the gravitational waves emitted at late time. When accretion dominates, 
the frequency and the amplitude of the gravitational- wave signal produced during the latest stages of 
the inspiral are nearly constant. In the exterior region we study a model in which the inspiral is driven 
by the emission of gravitational and scalar waves, which we describe completely at a relativistic level. 
Resonances in the energy flux appear whenever the orbital frequency matches the effective mass of 
the dark matter particle and they correspond to the excitation of the central object's quasinormal 
frequencies. Unexpectedly, these resonances can lead to large dephasing with respect to standard 
inspiral templates, to such an extent as to prevent detection with matched filtering techniques. We 
discuss some observational consequences of these effects for gravitational- wave detection. 
Subject headings: Dark matter — Gravitational-wave emission 



1. INTRODUCTION 

The Universe is populated with a plethora of different 
gravitationally-bound objects in dynamical equilibrium. 
Luminous, hydrogen- fueled stars like our Sun are sup- 
ported against collapse by radiation and gas pressure, 
whereas darker, compact and quiescent objects like neu- 
tron stars are prevented from full collapse by degener- 
acy pressure. Although evidence has been mounting for 
decades, only in recent years has it has become apparent 
that a completely different class of objects may, or must, 
also abound. Dark matter (DM) makes up a large frac- 
tion of galaxies and even though its exact nature is not 
known, it conglo merates into huge h alos around the cen- 
ter of galaxies (jBertone et al.lf2"005l ). Because all forms 
of matter gravitate, compact self-gravitating DM config- 
urations could therefore be a substantial component of 
our own galaxy. 

Compact DM objects have also been occasionally in- 
voked as an alternative to one of the most intriguing 
predictions of general relativity, namely the existence of 
black holes (BHs). Very massive main-sequence stars are 
dynamically unstable, a feature which is shared by most 
of the known compact configurations. Thus, "standard" 
stars cannot explain the dark, compact and supermassive 
objects lurking at the center of most gala xies, like the 
~ IO^Mq object in our own Milky Way (jGenzel et al.l 



Il996f ). It is widely accepted that BHs arc the most 
natural explanation for these supermassive compact ob- 
jects. Nevertheless, although actual observations sup- 
port the BH hypothesis, experiments showing the direct 
existence of an event horizon are still missing. In fact, 
some argue that an observational proof of the event hori- 
zon based only on elec tromagnetic observation is fun- 
damentally impossible (jAbramowicz "eFaII[200l . In an 
attempt to test the BH paradigm, exotic forms of mat- 
ter which possibly collapse to form supermassive hori- 
zonless objects, have been proposed. Besides their rel- 
evance for testing fundamental aspects of gravity, these 
objects may contribute to the da rk matter content of 
the Universe (jBertone et al.|[2005D . being thus relevant 
for particle physics and cosmology. The exact nature 
of the object at the center of our gal axy will soon be 
strongly constrained b y observations (iDoeleman et aP 
l200l:lFisi7eFaII [201II: faisenhauer et al.l I2008D . ' making 
this an exciting time to theoretically model and under- 
stand alternatives. 

Common approaches to probe DM in astrophysical set- 
tings are based on model-dependent interactions between 
the DM and the barionic sector. Such approaches usu- 
ally focus on the imprints these interactions leave on the 
evolution and equilibrium structures of astrophysical ob- 
jects. However, the equivalence principle guarantees that 
all forms of matter gravitate as predicted by Einstein's 
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general relativity, regardless of the (non)barionic nature 
of the constituent particles. As such, model-independent 
signatures of DM can arise from the study of gravita- 
tional self-interacting effects, like the existence of super- 
massivc DM configurations. See (jEda et al.ll2013f ) for a 
very recent proposal in this direction. 

It is widely believed that measurements of gravita- 
tional waves (GWs) from inspiralling stellar-mass ob- 
jects into supermassive cor npact objec t s will map the 
entire spacetime geometry (jRvanl I1995L Il997| ) and will 
carry imprints of the nature of the central object. Thus, 
GW measurements are in principle able to test the ex- 
istence of compact DM objects and to discriminate be- 
tween betwen_^Hs and other type s of horizonless ob- 
jects (jKesden et al.l [20051: [Pmlil [20121 ). 

With the above as motivation, we study distinctive fea- 
tures of the extreme mass-ratio inspiral (EMRI) around 
supermassive scalar-field configurations. In this system, 
a small compact object with mass fip £ (1, 10)Mq spirals 
into a supermassive object with mass M G (10^, 1O^)M0. 
The late-time inspiral of this syst em is of interest to 
future space-based GW detectors (jAmaro-Seoane et al.l 
I2012al[bl ). The typical orbital period is (102,10") 
seconds and low-frequency GWs are emitted in the 
(lO"'', 10-2) Hz frequency band. In the EMRI hmit, 
the inspiral can last tens to hundreds of years in the 
detector band. In one year of observation, millions of 
radians are contained in the signal, encodmg rich infor- 
mation about the spacetime dynamics ()Rvanlll995f ). For 
this reas on EMRIs are ex ceptional probes of strong-field 
gravit y (iGair et aI1l2012f), modifications of general rela- 
tivity (Cardoso ct al."201l UPani et al.l[20TT l: lYunes et al.l 
[2012: Yagi et al. 2012) and of the nature of supermassive 
objects. 

Here we argue that the GW signal from an EMRI can 
be also used to probe the existence of exotic fields that 
constitute the dark content of the Universe. If the cen- 
tral object is a BH, the ins piral terminates w ith a merger 
and subsequent rin gdown (jBerti et al.l 120091 ). Instead, if 
the central object is formed by some compact DM con- 
figuration that interacts very weakly with barionic stars, 
the EMRI proceeds also in the interior of the object, 
contributing significant amounts of signal-to-noise ratio 
to the signal (jKesden et al.ll200"5l ). EMRIs are relatively 
clean systems that can be described with great accuracy 
within a perturbative approach. During most of the in- 
spiral the stellar-mass object can be considered as a test 
particle moving on a fixed background. The timescale 
for merger is much longer than the orbital period and 
the evolution can be described by an adiabatic approx- 
imation. At each instant, we consider that the particle 
follows a geodesic of the background spacetime and the 
secular evolution of the geodesic parameters be computed 
by solving the linearized Einstein's equations. In this 
way one finds the inspiralling orbit and the correspond- 
ing gravitational waveform. This procedure takes into 
account the main dissipative effects of the back-reaction. 
A more detailed analysis, which would also co nsider con- 
servative effects (jPoissonI 120041 : lBarack|[2009() is beyond 
the scope of this work. 

Among other models for s clf-gravitating exotic fields 
(e.g. axion stars (iKolb fc Tkach cv 1993), boson- fermion 
stars (jHenriques et al.lll990l ). etc.). boson stars (BSs) 



are particularly interesting because they arise as sim- 
ple solutions of the Einstein-Klein- Gordon equations, 
without requiring any exotic matter other than a mas- 
sive bosonic field (for rev i ews on the subject see (iJetzei 
T99I iSchunck fc Mielkd [20031 iLiebling fc Palenzuela 



2012( )). BSs are compact stars configurations that may 



be th ought of as a natural realization of Wheeler's 
geons (|Wheelerlll955[ ) for scalar fields. Unlike the orig- 
inal geons, BSs can admit stable configurations which 
share many features with central galactic objects, with- 
out having sing ularities nor horizons (Torres ct al. 2003 
iGuzmaJ 120061 ). being indistinguishable from BHs in 
certain regimes. Formation of BSs has been stud- 
ied extensively in the literature (IFrieman et al] | 1988| : 



Ferrell fc G lcisei^'1989^'Frieman ct al. lll989l:lGra sso'[l990l 
Liddlc 1990: Tkachcv 19911: lLiddle"fc Madscr 



Madsen 

19921: ISchunck fc Mi elkc 200 31). T he recent discovery of 
the Higgs boson " ijAad ct al.ii20Tl ) is of course a further 
motivation to study thi s type of solutions. 

BSs can be classified (ISchunck fc Mielkel[2003l) accord- 
ing to the scalar potential in the Klein-Gordon La- 
grangian. Depending upon the scalar self-interactions, 
the maximum mass of a BS spans the entire range from 
one to billions of solar masses. 

Rather than working on a case-by-case analysis, here 
we focus on generic features that can leave a characteris- 
tic imprint on the waveform. The inspiral can be divided 
into two different regimes: the motion in the exterior of 
the supermassive object and the motion in the interior. 
In Section [2] we discuss EMRIs inside generic DM con- 
figurations. We show that the inspiral is mostly driven 
by accretion of the scalar field in to the small compact 
object and by dynamical friction (jChandrasekharl [19431 : 
I0strikeilll999l ). rather than by GW emission. We include 
these effects in a Newtonian analysis and compute the 
signal emitted in GWs. The signal is markedly different 
from that arising during the merger into a supermassive 
BH and it also deviates substantially from the evolu- 
tion obtained when accretion and gravitational drag are 
neglected (Kc sden et al.l[2005l ) . While our results of Sec- 
tion[2]are fairly generic, to discuss the outer evolution we 
need to specify some relativistic model. This is done in 
Section [21 where we describe the EMRI around a spheri- 
cally symmetric BS. We show that the evolution is driven 
by the emission of gravitational and scalar waves, which 
we describe at fully relativistic level. We show that, dur- 
ing a quasi-circular evolution, the energy flux can be res- 
onantly excited. This leads to a large dephasing with 
important observation consequences for GW detection. 
We conclude in Section SI by discussing possible exten- 
sions of our approach. We use the signature (— , -|-, -|-) 
for the metric and in most of the paper we adopt natural 
units h = c = G = 1, unless otherwise stated. 



2. GW-SIGNATURES OF EMRIS INSIDE COMPACT DM 
CONFIGURATIONS 

We are interested in discussing simple models 
that can capture the salient features of the inspi- 
ral through a DM medium in the extreme mass- 
ratio limit. There is a variety of situations in 
which such inspiral can occur. The most obvious 
one is the existence of self-gravitating, compact DM 
objects. Another is the ins piral around Kerr BHs 
surrounded by bosonic clouds ([Arvanitaki fc Dubovskvl 
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20llt lYoshino fc Kodanial [2011: IPani et alJ [2012al: 
Witek et al.ll2012D . Ve~recentlv. the GW signatures 
of intermediate mass BHs with D M mini-halos has also 
been considered (|Eda et al.ll2013f ). 

Despite the muhitude of models one can conceive, two 
of the most important generic features of these configura- 
tions are: (i) they interact with standard barionic matter 
by purely gravitational effects and (ii) they are typically 
supported by self-interactions of a massive bosonic field, 
whose mass can range between 0(10-^°) eV (or smaller) 
and 0(1) TeV (or larger). In fact, as far as purely grav- 
itational effects are considered, many dynamical aspects 
do not even depend on the nature of the DM particles 
(e.g. spin, mass, coupling constants,...), but solely on 
their mass-energy distribution. 

Although compact DM configurations usually require 
relativistic effects to support their selfgravity, many fea- 
tures of the inspiral can be captured by a simple Newto- 
nian model. In this section, we consider compact New- 
tonian objects which are characterized by some nonvan- 
ishing density profile p{r) when r < R and zero outside. 
The simplest model is a constant-density staiQ supported 
by purely radial pressure, but more realistic nonconstant 
profiles are also possible. The crucial point here is to 
assume that these objects do not couple with standard 
matter in any way other than through gravitational in- 
teractions. In particular, the small orbiting perturber 
can penetrate the stellar surface and move inside the ob- 
ject. Our main interest here is to understand this motion 
and the corresponding GW emission. 

2.1. Accretion- and gravitational drag-driven inspiral 

While the motion in the exterior is driven by radiation- 
reaction only, when the particle penetrates the stellar 
surface other effects must be taken into account, namely 
the accretion of the nonbarionic mass onto the small com- 
pact object and the drag force due to the gravitational 
interaction of the orbiting perturber with its own wake. 
In the context of EMRIs, these effects have been taken 
into account to study the imprints of matter sur round- 
mg supermassive BHs (iBarausse fc Rezzollall2008D on the 
GW emission. We will show that in the interior of com- 
pact nonbarionic matter configurations, these processes 
are actually dominant and the inspiral is mostly driven 
by DM accretion rather than by GW emission. 

2.1.1. Accretion: ColUsionless versus Bondi-Hoyle 

Accretion of the scalar field onto the small compact 
object produces external forces that contribute to the 
secular evolution. As long as the accreted mass is much 
smaller than the total mass of the orbiting object, the 
assumption of quasi-stationary motion should provide a 
fairly accurate description. Let us start by some simple 
estimates, assuming the accretion cross section is roughly 
the geometrical one of the small compact object. For 
head-on collisions, the small compact object of mass fip 
traverses the entire diameter 2R of the star, accreting a 

^ Constant density solutions exist - both in Newtonian gravity 
and in general relativity — also when the fluid is anisotropic, a s 
in the case of selfgravitating scalar fields I IDev fc Gleiserl 12003) . 
However, we will show that possible anisotropies in the fluid are not 
important to model the inspiral, because accretion and dragging 
arc mostly insensitive to pressure in the extreme mass-ratio limit. 
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Fig. 1. — Level plots corresponding to the ratio -Rp/f = 
1,10~*,10~* in the (ctdm, "^dm) plane, where (tdm is the total 
self-interaction cross section for DM and is the mass of the 

DM particle. In the left-uppermost part of the parameter space 
the radius of the small compact object i?p is much larger than the 
mean free path I and accretion occurs at the Bondi-Hoyle rate 
In the lower-rightmost part of the diagram i?p <^ I and accretion 
is governed by Eq. ((Sjl. Straight lines refer to fj,p/M = 10 ~^ and 
ii = 2M, but they have a simple scaling with the inverse of the 
mass ratio and with the compactness of the central object. 

tube of length 2R and radius Rp. Therefore, 



M.. 



accreted 



2 i?2 



0.02^ < 1 

M 



(1) 



where wc used Rp sa /ip, R ^ 10 A/ and we have as- 
sumed constant density. Thus, during a single passage, 
the accreted matter has a negligible effect for head-on 
collisions. On the other hand, during the inspiral from 
the surface, the orbiting object can accrete much more as 
it sweeps through the equatorial plane, tearing a disk-gap 
of width /ip and area ^ R^ . In this case we get 



-^'^accrctcd 3 A/ N 



(2) 



for compact central objects. Thus, inspirals have to be 
carefully considered. 

In a more rigorous treatment, dynamical effects and 
the nature of the small compact object must be included. 
Accretion is described by 



(3) 



where p is the density of the DM configuration, v is 
the modulus of the velocity of the small object with re- 
spect to distant static observers and tr is the accretion 
cross section. The latter strongly depends on the phys- 
ical processes involved in the accretion and on the na- 
ture of the perturber. If the small compact object is 
a BH, whose radius Rp is much smaller than the mean 
free path € = (ctdm'^)"^ (o^dm and n being the DM self- 
interaction total cross section and the particle density, 
respect ively), then a n approximate formu la is g = nR^ff 
where (jShapiro fc Teukolskv..l983;.Giddings fc Manganol 
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[200l) 



R 



Rr, 



off 



(4) 



is the effective capture radius. Therefore, as — >■ 0, any 
particle that forms the supermassive object will eventu- 
ally fall into the small BH even when orbiting at large 
distance. Clearly, there exists a cutoff distance given by 
the radius R of the supermassive star and which corre- 
sponds to a minimum velocity, Wmin = Rp/R- In the 
small mass-ratio limit Vj^in ~ Hp/M <C 1. Therefore, for 
any v > Vtuin we get 



irpRl 



Rp < 



(5) 



If the small compact object is a neutron star, the effec- 
tive cross section is the minimum between the geomet- 
rical cross section and sum of the cross s ections upon 
the iiidividual nuclei of the star (see, e.g., (jGould et al.l 
ll990tlBertone fc Fairbah3l2008D ). Provided the scatter- 
ing cross section between DM particles and the stellar 
nucleons is larger than 10~''^^cm^, Eq. ([S]) is still a good 
approximation for the accretion rate at nonrelativistic 
velocities. If the small object is a white dwarf, the ge- 
ometrical cross section is typically larger than the indi- 
vidual scattering contributions, so that the effective cross 
section i s given by the sum of the in dividual nuclei cross 
sections (jBertone fc Fairbairnll2008l ) . In this case the ac- 
cretion rate will depend on the details of the microphysics 
and on the type of DM particles that make up the super- 
massive object. However, white dwarfs are not compact 
enough to sustain tidal forces and are likely to be tidally 
disrupted during the latest stage of the external inspiral. 
Thus, wc shall restrict out attention to accretion into rel- 
atively slow BHs and neutron stars only, both governed 
by Eq. ©. 

On the other hand, if the radius of the object is 
comparable to or larger than the mean free path, 
Rp 3> £, then accretion becomes a macroscopic pro- 
cess and cohesion forces and matter compressibility 
must be taken into account ([Shapiro fc Teukolskvlll983l : 
iGiddings fc Manganol 120081 ). When the small compact 
object is a BH, this ty pe of accretion is described by the 
Bondi-Hoyle fo rmula (Bondi||195^ iBondi fc Ho"vIill944l: 
iShapiro fc feukolskv,198a ) 
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where Cs is the speed of sound and A is a number of or- 
der unity which depends on the details of the fluid. In 
our numerical simulations we have assumed A = 1, but 
the results depend on A very mildly. In Fig. [TJ we show 
the straight lines in the DM cross section-mass ((Tdm~ 
todm) plane corresponding to Rp/i = 1, 10^"*, 10~^ for a 
mass ratio [ipjM = 1O^M0 and R = 2A/. In the upper- 
leftmost part of the parameter space Rp 3> t and accre- 
tion occurs at the Bondi-Hoyle rate ([6]). In the lower- 
rightmost part of the diagram Rp <^ £ and accretion is 
governed by Eq. ([S]). In the rest of this paper, we shall 
consider the two regimes separately. 

2.1.2. Drag force 



Another important effect is the gravitational drag 
which results in a dynamical friction force on the small 
comp act object traveling through the DM distribu- 
tion (IChandrasekhaH [19431 : lOstriked [T999h . The gravi- 
tational field of the small perturber is felt universally, 
including by the DM making up the compact configura- 
tion. Thus, a portion of this material is "dragged" along 
the inspiral, being tantamount to a net decelerating force 
acting on the perturber. In the EMRI limit, at the scale 
of the small compact object, the density of the medium 
is nearly constant, so that we can adopt the theory of 
dynamical friction for motion through constant-density 
media. For linear m otion, the dynamical force friction 
reads ([Ostrikeiill999| ) 
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9 ) 



with 



-v/Cs, V < Cs 
log (t^) , V>Cs 



(7) 



(8) 



where r^in is the effective linear size surrounding the 
small compact object. 

The gravitational drag force in the case of circular mo- 
tion has been derived in ([Kim fc Kimll2007|) . where it is 
shown that the curvature of the orbit will bend the wake 
at large distances from the perturber. However, it is 
found that the subsonic motion is remarkably similar to 
the linear case, whereas the supersonic motion deviates 
from the linear case when the size of the small perturber 
is not parametrically smaller than the the orbital radius. 
Within the extreme mass-ratio assumptions, the size of 
the perturber is always much smaller then the orbital ra- 
dius, i.e. rmin ^ f, so that we expect the linear-motion 
analysis to approximate well the generic orbital case. We 
have confirmed this expectation by checking the results 
obtained using Eq. ((7|) with those o btained by imple- 
menting the fitting formulae given in ([Kim fc Kiml[200l 
for the gravitational drag force in subsonic circular mo- 
tion. 

If the central object is very compact, the speed of 
sound is comparable to the speed of light, and, in the 
nonrelativistic limit we focus on, the motion is likely to 
be subsonic. On the other hand, supersonic motion is 
in principle allowed close to the surface and depending 
on the particular model at hand. Supersonic motion 
produc es a sharp enhancement of the d rag force when 
vRi Cs (|Ostrikeiill999i : iKim fc Kimll2007[ ). In particular, 
the supersonic drag force depends on and effects re- 
lated to the circular motion must be taken into account 
if '"min ^ ©(r/lO). Even when v Cs, in the extreme- 
mass ratio limit rmin r, the results obtained for the 
linear motio n are in remarkabl e good agreement with the 
exact ones l[Kim fc Kimll200"7f ). In our models we have 
assumed r^-^ = RcS ^ r. 

Dynamical friction might be very important during the 
inspiral. Indeed, since the force due to accretion reads 
Fa = fj'pV, we obtain 



IF, 



DF| 



3 ■ 



pcolV 



isioniGss 



V ^ Cs 



(9) 



5 



for the accretion rate ([5]) and 
I^dfI 

^Bondi 




(10) 



for the Bondi rate ©, respectively. In Eq. © we have 
considered Rp ~ and we have reintroduced the 

speed of Ught c for clarity. Thus, in the case of colh- 
sionless accretion, if the motion is supersonic but non- 
relativistic, <C w ^ c, dynamical friction dominates 
over accretion. In the subsonic regime, if ~ C(c), 
the drag force is subdominant in the nonrelativistic limit 
and we expect the late stages of the inner inspiral to be 
dominated by coUisionless accretion effects. Nonetheless, 
depending upon the ratio c^/c, there exists a crossover 
regime in which both effects are equally important. On 
the other hand, if accretion is governed by the Bondi 
rate, it is comparable to the gravitational drag force in 
the subsonic regime and it is negligible in the supersonic 
one. The results presented in the next section are ob- 
tained by including both the accretion and the drag force 
during the inspiral and considering both types of accre- 
tion separately. 

2.1.3. Gravitational radiation backreaction 

In addition to accretion and gravitational drag, the 
motion inside compact DM configurations is also driven 
by gravitational radiation-reaction, similarly to the usual 
inspiral around compact objects in vacuum. Neglecting 
radiation-reaction, the motion of the perturber is gov- 
erned by the gravitational force Fgir) = m{r)^p/r'^. For 
simplicity, in this paper we consider the following mass 
function and DM density profile: 



i{r)^M(^^y, p{r) 



aM 



(11) 



but our results extend straightforwardly to more realistic 
profiles. Note that the equation above reduce to the case 
of constant-density DM configuration when a = 3 and to 
the vacuum case (briefly discussed in Appendix when 
a = 0. 

Neglecting radiation-reaction, the energy and angu- 
lar momentum of th e perturber are conserved quanti- 
ties (Hda et all l2013[ ): 



1 



■2 



2/ipr2 (1 - a)i?"ri-" ' 



(12) 



L = fipr'^ip. (13) 
In the case of circular orbits of radius r, we get 



A/(l + a) 
''^i?"ri-"2(a-l) ' 



and the Keplerian frequency reads 



m{r) 



(14) 
(15) 

(16) 



For circular orbits the radius and the orbital frequency 
are constant and, therefore, these orbits dissipate energy 



throu gh the standard quadrupolar formula (jMaggiord 
[200l : 



GW 



(17) 



where we have used Eq. ([16)) . Finally, differentiating 
Eq. ([H]) and using the energy balance law, E ~ — i?GWi 
we obtain the evolution equation of the orbital radius 

64MVp 



5(a l)r3-2"i?2a ' 

whose solution reads 

r{t) = R{l+t/Ta)~^i^ 



(18) 



(19) 

with (a - 2)to = 5(1 + a)R'^ /{USAP^ip). Note that the 
late time behavior of the solution above strongly depends 
on a: for a < 2 the radius drops to zero in a finite 
time I To I, whereas for a > 2 the radius approaches zero 
asymptotically. In the singular case a = 2 the behavior 
is exponentially suppressed. 

Let us now show that radiation-reaction effects are 
small compared to accretion. Simply by angular mo- 
mentum conservation and using Eq. (|16|). we obtain, for 
the secular evolution of the radius under accretion, 



"(t)/ip 



1 + a fip (r) 

By comparing with Eq. (jlSp . we obtain 

s 3/9 . . 3(o-l) 

^•collisionlcss _ 5q; / _R y / R 
rreaction 32 \Mj \r{t) ^ 

. ^. 1 

^Bondi 5aA R { R 



^reaction 

32c3 M \r{t) 



(20) 



(21) 



(22) 



for coUisionless and subsonic Bondi accretion, respec- 
tively. Interestingly, in both cases and for any a > 1, 
the late-time inspiral in the interior of DM compact con- 
figurations is generically dominated by accretion and not 
by radiative dissipation. 

2.2. Numerical evolution of the inspiral in the interior 

In this section we describe the secular evolution of the 
small perturber inside a spherically-symmetric DM con- 
figuration as driven by DM accretion and gravitational 
drag. As proved in the previous section, GW radiation- 
reaction is a small effect compared to accretion, so we 
can safely neglect it here. In Newtonian theory, the 
accretion- and gravitational drag-driven inspiral is de- 
scribed by the following system: 



(j-pr + lJ^p{r 
rfiptp + ^J.p{r(f + 2r(p) 



FBF,r , (23) 

(24) 



together with Eq. or Eq. (jB]) for the evolution of fJ,p{t). 
In the equations above i^DF.r and F-Q-p,^ are the two com- 
ponents of the gravitational drag force. Namely, 

r rip 

^DF,r = -Fbp- , ^DF,i^ = ^DF ■ (25) 

V V 

where v'^ = + r'^ip'^ and Fuf is the linear dynamical 
force friction, Eq. As discussed above, using the 
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formula for linea r motion is well jus tified in the extreme- 
mass ratio Hmit (|Kim h Kimll2007[) . 
Equation p4p can be directly integrated in two ex- 



tremal limits. Neglecting gravitational drag, and for any 
accretion rate, we get 



(26) 



Fig. 2. — Secular evolution of the orbital parameters of a point particle orbiting a constant density, Newtonian star with radius R = 6M. 
The particle starts at r(0) = 8M with initial eccentricity = e{t = 0) = 0.1. When r > R (blue curves), the evolution is radiation-driven 
through the quadrupolc formula (cf. Appendix [A]l . When r < R (red curves) the evolution is driven by dynamical friction and by (i) 
collisionless accretion (left panels, cf. Eq. l(5j) or ii) Bondi-Hoyle accretion (right panels, cf. Eq. lO). Upper panels: Cs = 0.6 and the 
motion is always subsonic. In the interior the orbits circularize quickly. Lower panels: Cs = 0.2; the inspiral in the interior starts supersonic 
and the evolution is dominated by dynamical friction. When v < Cs, the evolution proceeds qualitatively as in the upper panel. In the 
small left panels we show (from top to bottom): the radial position in polar coordinates, the module of the particle's velocity and the time 
evolution of the mass-ratio. The evolution does not qualitatively depend on the accretion rate formula used. 

the parameters and for nonconstant density profiletQ. 

In absence of dissipative effects, the motion inside 
constant-density distributions is worked out in Ap- 
pendix |b1 where we show that the small body moves 
on ellipses centered at the origin. The quasi-elliptical 
evolution is then governed by Eqs. Figured] 
summarizes the quasi-elliptical evolution of the orbital 
parameters in the subsonic (upper panels) and in the su- 
personic (lower panels) cases, both for collisionless accre- 
tion (left panels) and for Bondi accretion (right panels), 
respectively. In the top panels, we show the evolution of 
the orbital parameters for the case Cg = 0.6 (in units of 
the speed of light), i.e. the motion in the interior of the 
object is always subsonic. In the large panels we show 
the orbit in cartesian coordinates. The small panels refer 
to the evolution of the orbital radius, the velocity and the 
mass of the small object, respectively. To help visualiza- 
tion, we have considered the unrealistic initial mass ratio 
^p{t = 0) = 10~^M and the number of cycles grows with 
the inverse of fip/M. In the interior, accretion is very ef- 
ficient and the mass-ratio becomes of order unity during 
the inspiral, so that our approximation breaks down. 
In the interior (red curves) the orbit circularizes, re- 



L 



which can be also obtained from the conservation of the 
angular momentum L. On the other hand, if we neglect 
accretion (/ip =const) and in the limit ?j <C c,,, we get 



rip 



L 



■ cxp 



Airp 
3Z 



^ipt 



(27) 



where we have assumed constant density. 

In the general case, the system has to be 

integrated numerically for a given density profile p{r) 
and some initial conditions. The latter are chosen at the 
time the particle reaches the radius of the object at the 
end of the quasi-elliptical, radiation-driven inspiral in the 
exterior. 

Some results arc shown in Fig. [5J As an example, we 
have considered a constant-density, Newtonian star with 
radius R = 6M and a point-like particle located at r = 
8M at t = 0. The initial eccentricity is = e{t = 0) = 
0.1. Similar results can be obtained for other choices of 



^ Nonetheless, in the deep interior of stellar configurations the 
density is nearly constant, so that we expect constant-density pro- 
files to be a good approximation for the latest stages of the inner 
inspiral inside more complicate matter configurations. 
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gardless the details of the accretion rate. This can be 
proved analytically in the case of accretion-driven inspi- 
ral. Indeed, using the results of Appendix |B1 together 
with conservation angular momentum conservation and 
the energy balance law, it is easy to derive the secular 
evolution of the semi-major axis b and of the eccentricity 
in the small-eccentricity limit: 



^accretion ^ 



^accretion ' 



4Mp 



(28) 
(29) 



The equation above shows the important result that cir- 
cular motion remains circular. Because /ip > 0, or- 
bits which are slightly eccentric will tend to circularize. 
Therefore, as clearly shown in Fig. [31 circular inspiral is 
an attractor of the motion. Remarkably, this result is 
independent of the specific accretion process. 

In the bottom panels of Figure [H we show the same 
quantities as in the top panels but for the case Cg = 0.2, 
i.e. the motion in the interior starts supersonic and, as 
the velocity decreases, the inspiral enters the subsonic 
regime. As expected, during the supersonic phase the 
motion is dominated by dynamical friction and the ve- 
locity decreases abruptly. As the particle enters the sub- 
sonic regime, the evolution proceeds qualitatively as in 
Fig. [2l However, the inspiral in the supersonic case is 
much faster and the orbits do not have time to circu- 
larize. Note that the evolution does not qualitatively 
depend on the accretion rate formula used: the results 
shown in the left and right panels of Figure [3] are quali- 
tatively similar. 

Finally, note that the accretion can be extremely ef- 
ficient. In the latest stages of the inspiral, the orbiting 
object will start accreting very fast [cf. (|36|) below] and 
an arbitrarily small initial mass will accrete an amount 
of matter M in a finite amount of time. Therefore, our 
small mass-ratio hypothesis will eventually break down. 
If the orbiting object is a neutron star, during the inspi- 
ral it will devel op a DM core and it will l ikely collapse 
to form a BH (|Bertone fc FairbairnI 120081 ). The signal 
emitted during the collapse is another distinctive feature 
that must be quantified by a relativistic analysis. 

2.3. Gravitational waveforms 

Here we present the waveforms emitted during the in- 
spiral shown in Fig. [21 Once the orbital radius, angular 
velocity and perturber mass are obtained as functions 
of time, we can use the standard quadrupole formula to 
compute h^(t) and hx (t) for a distant observer located at 
f wit h an angle view of (t, 0) [cf. Eq. (3.72) in (jMaggiord 
[200i]. 

In Figure [21 we consider l = tt/2, so that only h+ is 
nonvanishing and the wave is linearly polarized. The 
motion in the exterior is driven by the classical radiation- 
reaction mechanism, briefly summarized in Appendix 
for completeness. In panel a) we have neglected accre- 
tion and dynamical friction, so that the evolution pro- 
ceeds only through the radiation-backreaction. In panel 
b) we have considered accretion but neglected gravita- 
tional drag and the signal has a constant amplitude and 
constant frequency. We discuss this case in detail in the 



next section. In panels c)-c) we have included both accre- 
tion and dynamical friction with different constant values 

of Cs- 

The last panel shows the supersonic case described in 
the lower set of panels of Figure [2l As expected, the 
contribution of dynamical friction becomes dominant as 
Cs ^ c and we observe two effects: (i) the amplitude 
of the signal decreases in time and (ii) the total time of 
the inspiral in the interior quickly decreases for smaller 
values of c^. The latter effect is simply due to the extra 
dissipative channel during the evolution. 

These features would be completely missed by an evo- 
lution purely driven by radiative effects. Our analysis 
does not include relativistic effects, but provides a strong 
case for including accretion and gravitational drag effects 
in the inspiral inside compact DM configurations. 

2.4. Analytical Fourier waveforms in the 
stationary-phase approximation 

Some of the qualitative features of the waveforms pre- 
sented in Figure [21 can be analytically captured by a toy 
model in which we consider a small object of mass Hp on 
a quasi-circular orbit inside a Newtonian star and whose 
secular evolution is driven by a single dissipative effect. 
At Newtonian order, it is possible to obtain analytical 
templates of the waveforms in Fourier domain in the case 
of accretion-driven and radiation-driven inspiral. To low- 
est order, /ip, r and Q are constant and the Newtonian 
waveforms simply read: 



'1 



- COS^ L 



hx[t)^ cos L. sm(LJGwt) 

c^r 



cos{ujGwt) , (30) 
(31) 



where wqw = 2f2. Then, dissipative effects can be in- 
cluded by replacing the constant parameters wqw: ^ and 
/ip by a;Gw(i)i i^it) and fip{t), where the secular time 
evolu tion is governed by the specific dissipative mecha- 
nism (|Maggiorel2008f ) . In the next sections, we shall treat 
collisionless accretion, Bondi accretion and radiation- 
reaction separately. 

2.4.1. Collisionless accretion, Rp <^ i 

In order to isolate the effects of accretion, let us neglect 
GW reaction and gravitational drag. As we discussed 
above, if Cg = 0{c), the gravitational drag is a small 
effect compared to accretion, at least in the nonrelativis- 
tic regime and for collisionless accretion. Let us then 
consider the orbital evolution driven by accretion only. 
Because the system is in isolation, the total angular mo- 
mentum is constant, even when accretion is included. On 
the other hand, the binding energy evolves during accre- 
tion. A simple and powerful result can be obtained solely 
by angular momentum conservation. 



l^p{t)ritf 

For quasi-circular orbits, v 
quency reads as in Eq. (|16p . 



L 



(32) 



= rfl and the orbital fre- 
Using Eq. (|32p . we obtain 



L 



sjGm[r)r[t) 



(33) 
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Fig. 3. — Gravitational-wave amplitude fc_|_f along (t, 0) = (7r/2, 0) as a function of time when collisionlcss aceretion (left panels) and 
Bondi accretion (right panels) arc considered, respectively. Small panels refer to: a) usual radiation-driven inspiral (neglecting accretion 
and dynamical friction); b) accretion-driven inspiral (neglecting radiation and dynamical friction); c) inspiral driven by accretion and 
dynamical friction (neglecting radiation) with Cs = 0.6, which corresponds to the orbits shown in the top panels of Fig. [2] d) same as panel 
c) but with Cs = 0.5; c) supersonic regime with Cs = 0-2, which corresponds to the orbits shown in the bottom panels of Fig. [2] Remaining 
parameters as in Fig. [2] 



where, here and in the rest of the section, we restore 
factors G and c for clarity. FinaUy, using Eq. ([5]) we get: 



St:LG 



p{r)r{t) 



m{r) + ATip{r)r{tY 



Rp<t:e. 



(34) 



Once a density profile is specified, the equation above can 
be solved for r{t). Then, the other dynamical quantities 
fi{t) and fl{t) respectively read as in Eqs. ((55)) and 

To be concrete, let us consider the density profile given 
in Eq. (|TT|) . In this case, Eq. (p4|) can be directly inte- 
grated: 



r(t)=i?(l-tAinsp)^ 



(35) 
(36) 
(37) 



where p.p is the mass of the particle at the time t = 0, 
with R = r(0) and il; ~ ^l(O), and we have introduced 
the duration of the inspiral, 



6GaL 



87rG3/2a(p)//^*) 



where (p) = MI / {AttR^). 

With this solution at hand, we c an compute th e 
corresponding gravitational waveforms ()Maggiorell2008[ ) . 
Standard treatment allows to write Eqs. ([30)) and (|3T)) to 
lowest order as 



h+ {t) = (tret ) cos W{t^at ) , 
hx{t)=Ax (tret) sin w{tret) , 



(39) 
(40) 



where tret is the retarded time and w 
the case at hand, we get 



12ti, 



GM 



i?3 \t 



2jl dt'n{t'). In 



(41) 



'insp 



where r — tinsp ~ t and rog — w{tinsp)- Note that the 
orbital frequency 51 is constant when the density is ho- 
mogeneous. When the density is nonconstant it is rele- 
vant to compute the Fourier transform of the waveform. 
Adopting a stationary phase approximation (jMaggiord 
[2008h . we obtain 



where 



A+e 



2*4 



Axe 



(42) 



f-insp 



1 



c4f y 2(a - 3) 



Ax — 



2 cos L 
1 + cos^ 

1 It- 

3 + a 



■A 



+ 1 



R^ V 
GMJ 



R ^ I ^ 



-2^/ [tc + 
11/2. 



IT 

4 



(43) 
(44) 



(45) 
(46) 



where / = 2ri/(27r). Therefore, for a > the amplitude 
and the phase increase with the frequency. 
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2.4.2. Bondi-Hoyle accretion, Rp 3> 
In this case, from Eq. ^ we get 



SttLGX p{r)y^Gm{r)r{r) 
m(r) + 4Trp{r)r{ty 



Rp->1, (47) 



which, for the density profile given in Eq. (jlip . can be 
directly integrated 



insp 



(48) 
(49) 
(50) 



where we have assumed a ^ 7 and 
jB) _ 4MRy\l + a) _ 



insp 



(7 - a)aXL{GM)y^ AT:G^a{l - a){p)\y!^^ 



Similarly to the previous case, using the solution above 
we can perform a stationary phase approximation to 
compute the amplitude and the phase in the frequency 
domain. The relevant ones read 



2i 



G- 



l + g 3q-1.3 15-o 



r 



l/("-3) 



2 insp 



d3 \ 2(a-3) 

gm) 



~27r/ ( + - 



(51) 



(52) 



2.4.3. Radiation-reaction 



The motion driven solely by GW radiation reaction has 
been worked out in Section 12.1.31 Using Eqs. ([T^ and 
IH]) we obtain, to lowest order, 



A, 



(0 



To (a — 2) 1 + cos^ L 
a-3 2 



G' 



l/(a-3) 



5 — 3a 



2) fGM\ 



■2(o,-2) 
a-3 



Wo 



TT 

4 



(53) 



(54) 



where tq is defined below Eq. ([T! 

Comparing the waveform amplitudes and phase in 
Fourier space shows that each specific dissipative mech- 
anism produces a very peculiar signal. In particular, the 
signal produced by accretion-driven inspiral is dramati- 
cally different from that arising in the classical radiation- 
driven inspiral, already at Newtonian level. Our analysis 
can be easily extended to more realistic density profiles. 
In general, the signal strongly depends on the DM pro- 
file and the waveform parameters are also sensitive to 
the type of accretion process. Thus, the GW signal can 



be used to constrain the nature and the properties of 
compact DM configurations. 

3. RELATIVISTIC ANALYSIS FOR THE EXTERNAL 
INSPIRAL 

A self-consistent relativistic analysis of the inspiral 
around scalar-field compact configuration is fairly in- 
volved already at first order in the mass ratio. Indeed, 
due to the coupling between gravitational and scalar per- 
turbations and to the presence of a background scalar 
field, a small perturber will also source scalar waves, even 
if it is formed by purely barionic mass. On the other 
hand, this property allows for a much richer phenomenol- 
ogy that would be missed if a fully relativistic analysis 
is not properly performed. To be concrete, in this sec- 
tion we focus on specific models for selfgravitating DM 
objects. We investigate the emission of gravitational and 
scalar waves sourced by a test-particle in circular orbit 
around a boson star (BS). In this section we discuss the 
main features of this process. The mathematical pro- 
cedure is standard but technic ally involved. A det ailed 
analysis will appear elsewhere (jMacedo et al.|[20T3|) . 

3.1. Relativistic models of supermassive DM objects 

For concr eteness, in th e following we focus on BS con- 
figurations (jLiebling fc P alenzuela 2012]) which are rel- 
ativistic solutions of the Einstein-Klein- Gordon theory 
but most of our results hold at a qualitative level for 
different models. 

We consider the Einstein-Klein-Gordon theory: 



S 



R_ 

2k 



-.g''''aa$*96$-14(|$|' 



where k = Stt and S'mattor denotes the action of any bar- 
ionic matter field. The Einstein-Klcin-Gordon equations 
read 



-Rah 



1 



-gabR' 



T. 



ah 



Tnmattcr'i 
^ab ) 



1 



--da (V^S^'S,*) 



(55) 
(56) 



V -3 
where 

= da^*dk^ + 4**9a<i> - 9ab (9^**9^$ + Kd^H) , 

is the energy-momentum of the scalar field. For a com- 
plex scalar field, Eq. (|56p is supplied by its complex con- 
jugate. We will focus on spherically symmetric selfgrav- 
itating objects, whose line element is 



dsn 



_e''('-)di2 -f e"('')dr2 + r^{dB'^ + sin^ Odip^) , (57) 



whereas the background scalar field reads 

$o(i,r) = 0o(r)e-'"*, (58) 

with 4n){r) being a real function. The field equa- 
tions for u(r), u[r) and 4)o{r) form an eigenvalue prob- 
lem for u and they can be solved by standard mcth- 
ods. see e.g. (iKauol [19 68: Ruffini & Bonazzola 196£ 
Golpi et al.lll986t lGleise'i1il988i: ,Gleiser fc Watkinsi,198" 



Schunck fc Mielkell2003[ ). The equations can be recast in 
the form of Einstein gravity coupled to an anisotropic, 
nonbarotropic fluid, whose energy density and pressure 



10 







0o(r = O) 








Mu] 


R/M 


mini-BS I 
mini-BS II 




0.0541 
0.1157 


0.853087 
0.773453 


0.63300 
0.53421 


7.86149 
4.52825 


0.54000 
0.41319 


12.4194 
9.03368 


massivc-BS I 
massive-BS II 




0.0188 
0.0309 


0.82629992558783 
0.79545061700675 


2.25721 
1.92839 


15.6565 
11.3739 


1.86513 
1.53394 


6.9362 
5.8981 


solitonic BS I 
solitonic BS II 




0.0371 
0.0389 


0.1220326382068426831501347644 
0.1236174876926880171453994576 


7.36961 
6.77654 


22.8587 
20.2924 


0.89933 
0.83770 


3.1017 
2.9945 



TABLE 1 

BOSON-STAR MODELS USED IN THIS WORK. THE QUANTITIES M, R AND REPRESENT THE MASS OF THE SOLUTION, ITS RADIUS AND THE 
MASS OF THE SCALAR FIELD. FOR MASSIVE-BS AND SOLITONIC-BS CONFIGURATIONS WE USED THE FIDUCIAL VALUES A = SOOnlJ? AND 

(To = 0.05. The labels "I" and "11" respectively refer to stable and unstable configurations with respect to radial 

PERTURBATIONS. THE SIGNIFICANT DIGITS OF THE EIGENVALUE UJ DO NOT REPRESENT THE NUMERICAL PRECISION, BUT THEY SHOW THE 

FINE TUNING NEEDED TO OBTAIN THE SOLUTION. 



are defined in terms of the stress-energy tensor of the 
scalar field, T^. Namely, 

P^-T^t' = u;'e-^cj^l + e-{^',f + , (59) 



Pr — 1 r ^ LO e 

Pt = T e 



'"-"■^'j + e-^l^Jjf-K", (60) 



(61) 



where = Vs{(j)o) ^^'^ Pj Pr £md pt are the density, 
radial pressure and tangential pressure of the fluid, re- 
spectively. Using a standard shooting method, we have 
constructed spherically symmetric compact BSs which 
are solutions of three different models, presented in Ta- 
ble H 

3.2. Geodesies around boson stars 

Stellar-size objects gravitating around supermassive 
BSs have a small back-reaction on the geometry and to 
first order in the object's mass move along geodesies of 
the BS background. Accordingly, GW emission by such 
binaries requires a knowledge of geodesic motion together 
with the consequent perturbative expansion of the grav- 
itational field. Many features of the gravitational radi- 
ation can be understood from the geodesic motion, in 
which we now focus. We will also focus exclusively on 
circular, geodesic motion. The rationale behind this is 
that it makes the calculations much simpler while re- 
taining the main features of the physics. Furthermore, it 
can be shown that generic ecce ntric o rbits get circular- 
ized by GW emission in vacuum ([Peter s 19(54), on a time 
scale that depends on the ma ss ratio. 

We follow the analysis by (|Chandrasekhaii[l985f ). the 
formalism for a ge neric background is presented in 
(|Gardoso et al.l[200g ). Following previous studies, we as- 
sume that the point-partic le is not directly coupled to 
the background scalar field (jKesden et al.|[2005l : iGuzmanI 
|2006: Tor res 2002,). The conserved energy E, the angu- 
lar momentum parameter per unit rest mass L, and the 
orbital frequency of circular geodesies read 



(62) 
(63) 
(64) 



with A{r) = e" as defined in Eq. (|57p . A summary of 
the geodesic quantities corresponding to the BS mod- 




t V 2rc 



els in Table [T] is presented in Fig. |4l Up to the in- 
nermost stable circular orbit of a Schwarzschild space- 
time, r = 6M, geodesic quantities are very close to 
their Schwarzschild counterpart with same total mass, 
as might be expected for very compact configurations. 
However, the geodesic structure at r < 6A/ can be 
very different (jKesden et al.|[2005l ). A striking difference 
is that stable, circular time like geodesies exist for BSs 
even well deep into the star (|Torre!jl2003 [G uzma ^ [20061 : 
iKesden et al.ll2005h . 

For less compact configurations - namely mini BSs and 
massive BSs - stable circular geodesies exist all the way 
down the center of the star. This is an important feature 
because, if radiative effects and external forces are small, 
the inspiral will proceed through a secular evolution of 
these orbits. Furthermore, the orbital frequency is 
roughly constant close to the origin. This also implies 
that the circular geodesies deep inside the BS are non- 
relativistic, as the velocity measured by static observers 
at infinity vanishes as the radius approaches zero. In this 
case, our previous Newtonian analysis should provide re- 
liable results. Therefore, we expect that even the circular 
orbits in the interior of the object will be accessible to a 
quasi-circular evolution that, as we previously discussed, 
is mainly driven by accretion and gravitational drag ef- 
fects. 

More compact configurations such as solitonic BSs may 
show truly relativistic effects. For these models R ~ 3M 
and an innermost stable circular orbit exists at r « 6M 
with M^isco ~ 0.06804. This is to be expected, as 
the background scalar field is exponentially suppressed 
and the spacctime is very close to Schwarzschild outside 
the solitonic-BS effective radius. Like in the case of a 
Schwarzschild BH, there exists an unstable light ring at 
roughly r„i « 3A/. The unexpected feature is the pres- 
ence of a second stable light ring at r/_, together with 
a family of stable timelikc circular geodesies all the way 
to the center of the star. This is clearly a relativistic 
feature, which is similar to the case of uniform density 
stars. The latter may also present two light-ring and sta- 
ble circular time-like orbits in their interior, depending 
on their compactness. The bottom panels of Fig. 0] de- 
picts a uniform density star with radius R = 3A/. In 
this case, the two light- rin gs degenerate in the star sur- 
face. It has been argued (jKesden et al.l 120051 ) that, for 
these very compact models, the orbiting particle plunges 
when it reaches the innermost stable circular orbit and, 
due to radiation effects, the eccentricity of the orbits in 
the interior of the BS will increase. However, as we have 
discuss, accretion effects are dominant in the interior of 
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mini BS massive BS solitonic BS 




05 10 5 10 12345 



r/M r/M r/M 

Fig. 4. — Circular geodesic motion for different BS models and configurations (cf. Tabic \Tl ■ In the top, middle and lower row we show 
the orbital frequency f2, the energy Ec and the specific angular momentum Lc, respectively. Each column refers to a different BS model. 
From left to right; mini-BS, massive-BS and solitonic BS. For each model, we compare the geodesic quantities to those of a Schwarzschild 
BFI and for the solitonic BS model we also compare to the metric elements of a uniform density star with R = 3M. In the last column, the 
markers indicate the outer last stable orbit for solitonic BS configurations, which is approximately given by r fa 6Af and MQisco ~ 0.06804. 
The light-rings are given by rj^ RJ 2.72093A/ and r™, « 2.9812M, with A/Qj^ » 0.188818 and MUi^ » 0.192453, for the first configuration 
and ri_ RJ 1.91163Af and rm » 2.99883Af, with MQi_ as 0.184590 and MQi^ » 0.192452, for the second one. 

the star and they also contribute to circularize the orbit. 
Our results are based on a Newtonian analysis, which is 
nonetheless accurate well deep inside the star, where the 
velocity is small. It is therefore possible that, after the 
initial plunge, the particle will have access to these circu- 
lar geodesies, whose evolution is governed by accretion, 
rather than by radiative effects. 

3.3. Point-particle orbiting a boson star, resonant fluxes 
and quasinormal modes 

The gravitational and scalar energy fluxes emitted dur- 
ing the quasi-circular inspiral of a test-particle around 
a compact BS can be derived at fully relativistic level. 
The emission is governed by an inhomogeneous system 
of equations, whose regular solutions can be constructed 
via standard Green's function techniques. We present 
here the main res ults, whereas a det ailed analysis will 
appear elsewhere ()Macedo et al.ll2013D . 

The axial sector is governed by a single equation which 
does not involve scalar perturbations. Due to the explicit 
form of the source term, the axial flux is vanishing for 
even values oi I + m. In Fig. [SJ we show the dominant 
/ = 2, TO = 1 contribution of the axial flux for various BS 
models and compared to that of a Schwarzschild BH. 

The polar sector is described by an inhomogeneous sys- 
tem of coupled equations. A gener al method to solve this 
class of problems was pr esented in (jPani et al.ll2011| ) (see 
also () Molina et al.|[2010l) ). The contribution of the polar 
radiation to the total energy flux is nonvanishing only 
when I + m is even and it is maximum when I = m. The 
dominant contributions is shown in Fig. [6] as a function 




mini BS 

massive BS 

solitonic BS 

BH 



I I 



8 



Fig. 5. — Dominant ( = 2, m = 1 contribution to the axial 
gravitational flux emitted by a point-particle orbiting a BS for the 
stable BS configurations used in this work, compared to that of a 
Schwarzschild BH. The most compact configurations are closer to 
the BH case, and both solitonic configurations for r > 3Af have 
basically the same values of the BH case. 

of the orbital distance for mini BSs and massive BSs. 
The axial and polar sector present similar features: at 
large distance the deviations from the BH case are basi- 
cally indistinguishable, whereas for stable circular orbits 
inside the BS the energy flux quickly decreases as the 
orbit shrinks. In this limit, the orbital velocity is nonrel- 
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ativistic and our results agree with a simple quadrupole 
formula (|Maggiorell2008l ). 

However, the polar flux shows some sharp peaks 
which correspond to specific resonant frequencies. This 
interesting phenomenon occurs quite generically for 
smal l objects orbiting relativistic compa ct stars (see, 
e.g. (|Pons et al.ll2002l: iGualtieri et al.ll200H) l and it is re- 
lated to the excitation of the quasinormal mode (QNM) 
frequencies. For a point-particle with orbital frequency 
Q, the resonance condition reads 



mfi = an, 



(65) 



where m is the azimuthal number and an is the real part 
of the QNM frequency. In other words, if the characteris- 
tic frequency of the BS matches (multiples of) the orbital 
frequency of the particle, sharp peaks appear in the emit- 
ted flux. This can be modeled in terms of a simple har- 
monic oscillator, where the orbiting particle acts as an ex- 
ternal force and where the an is the proper frequency of 
the system. In this picture, the imaginary part of the fre- 
quency (7/ is related to the damping of the oscillator and 
it is roughly proportional to the width of the resonance, 
while the quality factor an/ai is p roportional to th e 
square root of the resonance height (iPons "eFaI1[200l . 
In agreement with this model, the resonances shown in 
Fig.inicorrespond to the lowest damped QNMs, which are 
presented in Table [2] The QNM spectrum of a BS is rich 
and comprises several classes of modes (|Yoshida et al.l 
[1991 . Here we have found a novel class of lowest damped 
modes, which can be excited during the inspiral due to 
their low frequency. A complete analysi s of the QNM 
spect rum will appear in a separate work (jMacedo et al.l 
[2011 . 



Model 




-ajM 




<5</>GW [rads] 


mini-BS I 


0.0757 


3 X 10-5 


0.01 


6 X 10=* 


massive-BS I 
massive-BS I 
massive-BS I 


0.0909 
0.1616 
0.2136 


6 X 10-^ 
5 X 10-S 
1 X 10-^ 


0.13 
0.02 
0.007 


9 X 10" 
2 X 10^ 
4 X 10^1 



TABLE 2 

Polar quasi-bound modes a = an + iai of mini-BS and 

MASSIVE-BS configurations CORRESPONDING TO THE 
RESONANCES SHOWN IN FiG. [6|f0R I = 2 AND COMPU TED BY A 
DIRECT INTEGRATION METHOD I|MACED0 ET AL.|[2013I ). We ALSO 
PRESENT THE HEIGHT OF THE RESONANCE NORMALIZED BY THE 
MASS RATIO, h/v'^ , AND THE TOTAL GW DEPHASING (S0GW AS 
COMPUTED IN EQ. ll68t FOR Tobs = lYR AND M = IO^Mq. 



As shown in Fig. [51 the resonant frequencies may cor- 
respond to a stable circular orbit located outside the 
BS radius (as for the rightmost resonance in the right 
panel of Fig. [6|) or may correspond to stable circular or- 
bits inside the BS (as in the mini-BS case shown in the 
left panel of Fig. [6]). While reson ant circular orbits also 
occurs around perfec t-fluid stars (iPons et al.ll2002}) and 
other BH mimickers (jPani et "al]|2010[) . the existence of 
resonant geodesies inside the compact object is peculiar 
of BSs, due to the absence of a well-defined surface and 
due to the existence of stable circular orbits inside the 
star. Similar results as those shown in Fig. [6] may be de- 
rived for other choices of the parameters and for higher 
values of I. 

The existence of these resonances is intriguing, be- 
cause they appear to be a generic feature of com- 



pact objects supported solely by the self-gravity of a 
scalar field. For BSs, this novel class of modes corre- 
sponds to the scalar perturbations being localized close 
to the BS radius and decaying exponentially at infinity, 
while the gravitational perturbations propagate to in- 
finity as plane waves. These mod e s are usually named 
"quasi- bo und" st ates (|Dolanll2007t iRosa fc DolanI I20T2I: 
iPani et aI1l2012b[ ) and they are supported by the mass of 
the scalar field. In fact, any sufficiently compact object 
can support this class of modes in its interior. Constant 
density stars can support bound-state modes (i.e. modes 
with purely real f requency) for massiv e scalar perturba- 
tions with I > (jMacedo et al.l [20131 ). In the case of a 
BS, these modes acquire a small imaginary part which is 
related to the coupling between scalar and gravitational 
perturbations: even if the scalar flux is zero for bound- 
state modes, part of the energy carried by the scalar held 
can be converted into gravitational energy that is then 
dissipated at infinity through GWs. This also explains 
qualitatively why the imaginary part of these modes is 
small (i.e. why the resonances are generically narrow), 
because the dissipation mechanism is not efficient. 

The overall structure of the resonances is fairly rich 
and it depends on the values of /, m and on the specific 
models. However, in line with the case of a Schwarzschild 
BH, the modes have a hydrogenic-like spectrum. In that 
case, the location and width of the resonances can be 
computed analytically in t he small mass limit (Dctweileil 
[T98a ICardoso et al.|[20Tl[ ). For a Schwarzschild BH the 
quasi-bound modes a = aii + iaj read 



a I' 



41 



2, ,2 



2{n + l + l) 



{21 



1)! 



M{l + l + n)^i^+i)Y[\ + i\'T[l + i\ 



where n > is the overtone number. Therefore, as a 
approaches an there is a multitude of modes that can be 
excited and their separation in orbital frequency vanishes 
in the large I or large n limit. In the same limit the imag- 
inary part (and hence the width of the resonances) of the 
modes decreases very rapidly, as shown by the last equa- 
tion above. Our results are in qualitative agreement with 
this behavior, although in the BS models we considered 
the scalar field mass is not small, cf. Table [TJ Indeed, in 
the opposite regime when /xM <gc 1 the imaginary part is 
exponentially suppressed, at least when the background 
spacetime is a BH (Zo uros & Eardley 19791). We have 
found the same exponential behavior in the case of soli- 
tonic BSs (which have ^iM = 0(10), cf. Table [I]). This 
makes it extremely challenging to compute the polar flux 
for solitonic BSs and the corresponding flux would show 
extremely narrow resonances. This is the reason why we 
do not show the polar flux for solitonic BSs in Fig. ([6]) 
However, the effective mass for a complex scalar field 
also depends on uj a nd, correspondingly , the resonance 
condition is shifted (jMacedo et al.|[2013[ ). Our analysis 
generically shows that the resonant frequencies are of the 
order 

S2rcs = , (66) 

m 

and the width of the resonances decreases quickly for 
large values of m and for the overtones. It would be 
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mini BS 



massive BS 




r /M 
p 



Fig. 6. — Main multipolc contributions, I, m = 2, 3 and 4, for the mini and massivc-BS configurations. 



interesting to derive an analytical formula similar to the 
one above, for BSs. Remarkably, the resonant frequeneies 
only depend on the parameters of the theory - namely on 
the scalar field mass and on w - so that possible detection 
of resonant fluxes can be used to discriminate among 
different models and to tell a BS from a supermassive 
BH. 

Let us now estimate the prospects of observing such ef- 
fects. A generic framework to study the de tectability of 
narro w resonances has been developed by (jYunes et al.l 
I2012f) . While the orbiting body crosses the resonance, 
the emitted energy flux increases by orders of magnitude 
and the inspiral proceeds faster. Therefore, the main ob- 
servational consequence is a dephasing of the GW signal 
induced by the orbital acce leration at the resonance. Fol- 
lowing (jYunes et al.l 120121 ). we model the energy flux as 
a top-hat functiorEI: 



E = Eo + h{t)H [(54, - (i - ^rcs)'] 



(67) 



where Eq is the flux in absence of the resonance which 
is well approximated by the quadrupole formula at large 
distance, h{t) is a time-dependent height of the resonance 
occurring aX t = tj-cs, litres is its duration and H is the 
Heaviside function. In our case i5tros is much shorter than 
the orbital period, thus the leading-order for mula for the 
deph asing induced by the resonance reads (jYunes et al.l 
I2OII) 



5h MAQrcs To 



bs 



GWl 



16l^2 (MnrcsyO/^ M 



(68) 



where h = /i(tros): = fJ-p/M ^ 1 and Afii-cs is the width 
of the resonance in the frequency space. We estimate 
Afires ~ 2(77 and f^rcs ~ cr_R, whcrc an and cr/ are the real 
and the imaginary part of the QNM frequency. Finally, 



^ Note tfiat, as shown in Fig. |H] and consistently with a simple 
harmonic oscillator model (jPons et al. 2002), the resonant flux con- 
sists of a resonance and an ant ires onance. Both can be modelled 
by a top-hat function as in Eq. I I67I I , but they would have opposite 
sign. However, the height of the antiresonance is of the order of 
the Newtonian flux or less, whereas the height of the resonance is 
large by orders of magnitude. Therefore, we can safely neglect the 
contribution from the antiresonance and focus on a single top-hat 
function as in Eq. II67I I. 



we obtain 



^8.6 X lO^rads 



lO^Mm 



M 



10- 



aiAI 



10-5 



0.1 



ohs 

lyr 

10/3 



(69) 



where we have rescaled all quantities by typical values as 
obtained for the peak of the resonance and for the QNM 
frequencies (cf. Tabled]). Space-based observatories like 
eLISA/NGO will be sensitive to v ariations in the GW 
phase of the order of one radian (jAmaro-Seoane et al.l 
l2Q12al ib[). The estimate above shows that the dephasing 
induced by the flux resonance can easily be larger by or- 
ders of magnitude. In fact, such large dephasing implies 
that a matched-filtering search for EMRIs that uses gen- 
eral relativistic templates would likely miss the signal or 
detect it but extracting completely wrong physical pa- 
rameters. Depending on the parameters of the model, 
the height and the width of the resonance can be smaller 
and, for very narrow resonances, the dephasing will be 
negligible. Nonetheless, the bottom line of our analy- 
sis is that the quasi-circular inspiral of small compact 
objects around supermassive BSs will leave potentially 
detectable imprints that would be missed if the central 
object is a priori assumed to be a BH. 

Once the energy fluxes are computed, the evolution of 
an EMRI can be modeled using a n adiabatic approxi- 
mation and a Teukolsky evolution (jSundararaian et al.l 
[2007; Hughes 2000). This procedure will include dis- 
sipative effects to all Post-Newtonian orders, but it is 
valid only at first order in the mass ratio. While such 
relativi stic approach (together wit h including self-force 
effects (jPoissonll200l IBarackll200l ) is crucial for gener- 
ating accurate templates, no qualitatively new features 
will arise with respect to the case of inspiral around a 
massive BH. Our goal here is to point out new effects 
that can be used as a smoking gun for scalar-field con- 
figurations. One effect is the appearance of resonances 
discussed above. Another effect is the existence of sta- 
ble circular orbits inside the BS and the possibility that 
the small compact object proceeds all the way down to 
the origin, without plunging. This was discussed in Sec- 
tion [5] within a much simpler - although more general - 
Newtonian analysis. 

4. CONCLUSIONS AND OUTLOOK 
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According to general relativity, if self-interacting fun- 
damental scalar fields exist in nature they may collapse to 
form compact self-gravitating configurations whose mass 
ranges from one to billions of solar masses, depending 
upon the scalar potential. Future GW detectors will be 
sensitive to the signal emitted by neutron stars and solar- 
mass BHs orbiting supermassive objects like those pow- 
ering active galactic nuclei. In this paper, we have inves- 
tigated several distinctive features of the inspiral around 
supermassive scalar-field configurations in the extreme 
mass-ratio regime. Rather than working on a case-by- 
case analysis, we focused on generic features that can 
leave a characteristic imprint on the gravitational wave- 
form and which are fairly independent from the micro- 
physics of DM particles. 

Working in a Newtonian approximation, we have dis- 
cussed the inspiral in the interior of very generic DM con- 
figurations. If the small compact object interacts purely 
gravitationally with the scalar field, its motion will be 
described by quasi-elliptical orbits whose secular evolu- 
tion is driven by DM accretion and by gravitational drag. 
These effects dominate the inspiral and are responsible 
for a peculiar GW signal. If accretion dominates over dy- 
namical friction, the signal has a nearly constant ampli- 
tude and nearly constant frequency at late times. This is 
markedly different from the classical plunge which would 
occur if the central object were a BH and it might be 
used to discriminate the supermassive objects in galactic 
nuclei and to probe DM. We have shown that, already 
at Newtonian level, the waveforms in the Fourier space 
are strongly sensitive on the DM density profile. 

As an aside application, our results may also be rele- 
vant to study th e GW signal fro m the inspiral of putative 
primordial BHs (|Hawking|ll971| ) in the interior of neutron 
stars or to study the ins piral around Kerr BHs endowed 
with bosonic clou ds (lArvanitaki fc Dubovskvl 120111: 
Yoshino fc Kodamall2012tlPani et al.ll2012aHWitek et al.l 

Secondly, the motion of the small compact object in the 
exterior of the supermassive configuration is driven by 
the emission of gravitational and scalar waves, which are 
coupled to each other. Due to this coupling, a barionic 
test-particle in quasi-circular motion can resonantly ex- 
cite the scalar QNMs of the central DM object. We have 
demonstrated this by considering some specific models 
of relativistic BSs in spherical symmetry. These reso- 
nances appear as sharp peaks in the energy flux emitted 
in GWs and would result in a faster inspiral, leading to 
a dephasing which can have observational consequences 
for future detectors. The resonant frequencies are largely 
insensitive to the details of the system, and they mainly 
depend only on the mass of the scalar particle. We have 
shown that the resonances correspond to the excitation 
of a novel class of QNMs, which have a much longer life- 
time and should therefore dominate the late-time signal 
during the gravitational collapse and during the ring- 
down. We have discussed the QNM spectrum and the 
linear response of the system at fully relativistic level. 

The plethora of DM candidates, modified gravitational 
theories and models for BH mimickers makes it manda- 
tory to select generic features such as those described 
here. More rigorous case-by-case analysis might be per- 
formed if some of these signatures are eventually de- 



tected. Indeed, our work can be extended in several 
directions. A rigorous treatment of the accretion and 
drag effects in the interior beyond the Newtonian ap- 
proximation is lacking. Relativistic effects might be im- 
portant at the interface close to the radius, where ve- 
locities may be moderately relativistic and the motion 
could be s upersonic. Spe cial relativity effects were in- 
cluded in (jBaraussd I2007D and it would be interesting 
to extend our work by using those results and also by 
including Post-Newtonian corrections to the accretion- 
driven inspiral. On the other hand, the late-time in- 
spiral close to the origin is intrinsically nonrelativistic 
and we expect our results to be accurate in that regime. 
Furthermore, a simple extension of our analysis is to in- 
clude noncircular motion at fully relativistic level (see 
e.g. (jjvlartcl 200i), where the same extension has been 
performed for the inspiral around a BH). Our results 
also show that a relativ istic analysis as that performed 
by (|Kesden et al.l [20051 ) should be extended to include 
the effects of accretion and dynamical friction. Finally, 
we focused here on nonspinning solutions, while spin will 
certainly play a crucia l role. I n particular, possible super- 
radiance instabilities (jCardos o et al. 2008), which affect 
horizonless and spinning compact objects, can be con- 
sidered together with their possible imprint on the GW 
signal emitted during the inspiral. 

While this study was in its last stages, a related work 
appeared in the literature (|Eda et al.ll2Q13| ) ; while the as- 
sumptions of both works are different, and our framework 
more general, the broad conclusions are both optimistic: 
gravitational waves can be an efficient tool to study DM. 
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APPENDIX 

APPENDIX A: RADIATION-DRIVEN INSPIRAL IN THE EXTERIOR 

Wi thin a Newtonian approximation the inspiral in the exterior of the central object is simple and we follow (jMaggiord 
I2008D . In absence of dissipation, the motion is elliptical with the mass M located at the focus and the eccentricity 
e and the semi-major axis h are conserved quantities. At lowest order, the secular evolution of e[t) and b{t) driven 
by the GW emission can be modelled by the simple quadrupole formula. Clearly, this approach neglects any truly 
relativistic effect like the resonances discussed in the main text. On the other hand, the quadrupole formula provides 
a good approximation at large distances, when the orbital velocity is nonrelativistic. At first order in the mass ratio 
lip/M <C 1, the equations for the orbital evolution are 

15 %4(i_, 2)5/2 (,1+ 304^ j ' (Al) 

5 63 (1 - e^f^ V 24 96 y ' ^ ^ 

(l + ecos^)2 /m 

^= l-e2 V6^' ^^^^ 

where Lp is the angle in polar coordinates (r, (p) and the motion occurs in the = 7r/2 plane. The evolution is obtained 
solving the system above with some initial condition for e, h and Lp. The radius and angular momentum of the orbit 
read 

Finally, the GW signal produced during the evolution is governed by the functions and /ix, which can be written 
in terms of r(t), (p{t) and of the relative position of the distant observer (sec, Eqs. (PH)) and (^1]) in the main text). 
Note that the adiabatical evolution is valid provided radiation-reaction effects are small and the orbits evolve on 
timescales much longer than a typical orbital period. This is indeed the case in the EMRI limit, at least when the 
orbital separation is large. Relativistic corrections to the quadrupolar formula are indeed important in the final stages 
of the exterior inspiral, close to the radius of compact stars. 

APPENDIX B: MOTION OF PARTICLES IN A SPRING-LIKE POTENTIAL 

In this section we discuss the motion of a small particle inside a Newtonian, constant density star assuming the 
particle interacts only gravitationally with the star. This configuration should be a good approximation during the 
latest stages of the inspiral inside compact scalar configurations, where the density profile is nearly constant. The 
motion of a body inside a constant density medium is described by a gravitational potential of the form 

V{r) = ~^3^JiJ,+^^lJ,r\ (Bl) 

Here 7 = M/{2R^), /? = 3M/(2i?), where M,R are the mass and the radius of the constant density object and /ip is 
the test particle's mass. In polar coordinates, the motion of the body is described by 



lipr'^ = L, (B2) 



1 .2 L 



r^r^ + ^ + y(r)=E, (B3) 
where L, E are the conserved angular momentum and energy parameter, respectively. We can re-express the above as 

where = fJ,p{E + Piip)/L'^ . Changing to z = 1/r^ — l/fc, we get 



1 dry 2^lpE 1 2^lpV{r) 2 1 2/x2^r2 



dip 



Az^ + AB^ , (B5) 



where 



L2 

The solution to this equation is z = _B cos 2ip. Thus, we get 



+ kB COS 2(f) . (B7) 
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In cartesian coordinates, this can be expressed as the standard equation for an ellipse: 



3i_ r_ 

where a and b are the semi-minor and semi-major axis and 



1, 



(B8) 



(B9) 



1 + fcS ' " 1-kB 

Thus, we get the interesting result that test masses follow ellipses which are centered at the origin and whose eccentricity 
reads 

e= o?/b'^ = ^/2kB/{l + kB) . (BIO) 

Recalling that = b^{l ~ e^), we can re-express the energy and angular momentum in terms of the eccentricity and 
semi-major axis b as 



E = -Pfip+-/fipb\2-e''), 



(Bll) 
(B12) 



Finally, using A ~ ~ L/{2iip) and integrating over one orbit we get A ~ LT/{2^p), with T the orbital period 

and A the area of an ellipse. Using A — nab, we find the analog of Kepler's third law: 



2^2 



7(1 -e2) 



(B13) 



and we obtain the interesting result that the orbital period is completely independent from semi-major axis, as could 
be anticipated from the circular orbit case. Note that T cx p^^^^, where p is the density of the medium. 
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